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(57) ABSTRACT 

The present invention is a method for identifying unknown 
parameters in a system having a set of governing equations 
describing its behavior that cannot be put into regression 
form with the unknown parameters linearly represented. In 
this method, the vector of unknown parameters is segmented 
into a plurality of groups where each individual group of 
unknown parameters may be isolated linearly by manipula- 
tion of said equations. Multiple concurrent and independent 
recursive least squares identification of each said group run, 
treating other unknown parameters appearing in their regres- 
sion equation as if they were known perfectly, with said 
values provided by recursive least squares estimation from 
the other groups, thereby enabling the use of fast, compact, 
efficient linear algorithms to solve problems that would 
otherwise require nonlinear solution approaches. This inven- 
tion is presented with application to identification of mass 
and thruster properties for a thruster-controlled spacecraft. 

15 Claims, No Drawings 
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MULTIPLE CONCURRENT RECURSIVE 
LEAST SQUARES IDENTIFICATION WITH 
APPLICATION TO ON-LINE SPACECRAFT 
MASS-PROPERTY IDENTIFICATION 

5 

This application claims the benefit of the filing date of the 
previously filed provisional patent application, No. 60/468, 
411, filed on May 6 , 2003. 

FEDERALLY SPONSORED RESEARCH 10 

This invention was made with Government support under 
contract number NAS2-00065 awarded by NASA. The 
Government has certain rights in the invention. 

15 

COMPUTER PROGRAM LISTING COMPACT 
DISK APPENDIX 

Selected files used to implement a simulation of the 
invention, as applied to a thruster controlled spacecraft, are 20 
provided in MATLAB m-code on CD-R. All code is copy- 
right Edward Wilson. One original and one identical copy 
are provided. Machine format: IBM PC/XT/AT, or com- 
patibles. Operating system compatibility: MS-Windows. 
Line Terminator: ASCII Carriage return plus ASCII Line 2 s 
Feed. Control codes: none. Compression: uncompressed 
data. A printed listing of the files on the CD-R is also 
provided as an appendix to this specification. 

BACKGROUND 30 

1. Field of the Invention 

The present invention relates generally to information 
processing systems that are used to identify unknown 
parameters in a system, and more specifically to an infor- 35 
mation processing system that monitors motion-related sen- 
sors to accurately identify spacecraft mass-properties. The 
general field is commonly known as System Identification 
(ID). 

2. Prior Art 4 n 

The use of linear least squares regression for the identi- 
fication of unknown system parameters has been used and 
studied extensively, as by Lawson, C. and Hanson, R., in 
Solving Least Squares Problems, 1974, and Ljung, L., in 
System Identification, Theory for the User , 1999. However, 45 
the requirement that a regression equation be formed with 
the unknown parameters linearly represented, limits its 
direct applicability to many important problems, including 
the spacecraft application presented here. 

Ljung and several other authors have developed 50 
approaches for identification of parameters in highly non- 
linear systems using methods such as gradient-based opti- 
mization and neural networks. However, these methods are 
significantly more computationally intensive than those for 
linear problems. 55 

The remaining prior art items relate more specifically to 
the spacecraft mass-property ID aspect of the invention. In 
that problem, some of the unknown parameters multiply 
each other in the governing equations, resulting in a specific 
type of nonlinearity. 60 

Tanygin, S. and Williams, T., in “Mass property estima- 
tion using coasting maneuvers,” Journal of Guidance, Con- 
trol, and Dynamics, 1997, developed a least squares (LS) 
based algorithm to identify mass properties for a spinning 
vehicle during coasting maneuvers. The restriction to the 65 
case of a spinning spacecraft with no applied torques or 
thrusters firing limits its applicability considerably — either 
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to spacecraft that normally exist in this state, or by requiring 
other spacecraft to attain this state. 

Bergmann, E., et al., in “Mass property estimation for 
control of asymmetrical satellites,” Journal of Guidance, 
Control, and Dynamics, 1987, developed an ID approach 
using a Gaussian second-order filter as presented more 
generally by Gelb, A, et al., in Applied Optimal Estimation, 
1 974. The second order filter resembles an extended Kalman 
filter, but has extra terms to address the second order effects. 
This is significantly more complex and computationally 
intensive (by about two orders of magnitude) than the 
approach presented here, and may not produce better results 
for most spacecraft. The extra complexity may make it more 
susceptible to noise and parameter variations than the pre- 
sented methods. It assumes perfect knowledge of thruster 
properties. 

Wilson, E. and Rock, S. M., in “Reconfigurable control of 
a free-flying space robot using neural networks,” Proceed- 
ings of the American Control Conference, 1995, developed 
an ID method based on exponentially weighted RLS using 
accelerometer and angular rate sensors. The acceleration 
created by each thruster (reflecting both mass and thruster 
properties) was identified. This approach (identifying 
thruster acceleration rather than separately identifying mass 
and thruster properties) is more direct (since thruster accel- 
eration is the real value of interest from a control, estimation, 
or FDI standpoint), and probably better for vehicles with 
properties that are truly unknown (such as for the case where 
deflected thrusters are allowed, such as on the vehicle tested 
in that research). However, for most vehicles, certain prop- 
erties are well known, such as the thrust directions and 
locations in the structural frame. The present invention can 
take advantage of that knowledge to get better estimates of 
the properties that are not well known. 

BACKGROUND— OBJECTS AND 
ADVANTAGES 

In System ID, there is presently an extremely large step 
taken in terms of computational complexity when tackling a 
nonlinear vs. linear problem. The present invention 
addresses this for a broad class of problems, providing a 
novel approach to enable nonlinear system ID problems to 
be solved with linear methods. 

With regard to the spacecraft mass-property ID problem, 
due to the very small forces and torques present on space- 
craft on orbit and in free space, their mass properties are 
important from a control and estimation standpoint. Space- 
craft mass properties can only be calibrated with limited 
accuracy during ground testing, and change further once on 
orbit due to expulsion of fuel mass, reconfiguration (of 
antemiae, etc.), and for servicing robotic spacecraft, poten- 
tially variable payloads. Accurate ID of mass properties has 
been studied extensively as reported in the prior art, but 
those methods have not yet been implemented and tested on 
an actual spacecraft, possibly due to the computational 
complexity involved. 

The present invention was originally developed to address 
the spacecraft mass-property ID problem, and then extended 
to the general case. It enables highly accurate mass-property 
ID and is very simple, compact and fast. It has been 
implemented on an experimental spacecraft, tested on-board 
in zero-g aircraft testing, and is presently awaiting launch for 
space-based validation. 

Specifically, the primary object of the invention is an 
algorithm that enables the application of very powerful, fast, 
and compact system ID algorithms to a class of problems 
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that prior to this required solution with much more compu- 
tationally intensive, and fragile nonlinear approaches. It has 
been found to provide near-ideal results in testing on a 
realistic, important, and well-formed example application, 
where ideal indicates that the approximation error made in 
allowing this nonlinear problem to be solved with linear 
methods results in a negligible loss of accuracy in the ID 
result. 

BRIEF SUMMARY OF THE INVENTION 


4 

equation would be equal if noise were not present. The LS 
ID solution, x, minimizes the sum of the squares of the 
elements of the error, Ax-b. If the problem at hand can be put 
into regression form, with noise appearing only in the e term, 
5 x can be solved directly (i.e., this is a closed-form solution, 
rather than an iterative optimization as might be required if 
the equations can not be put into this form) using one of the 
following equations: 

Unweighted, batch algorithm: x=(A T A)~ l A T b (3) 


The present invention is a method for identifying 
unknown parameters in a system having a mathematical 
model describing its behavior and one or more measure- 
ments that are sampled regularly, and where the system’s 
governing equations cannot be manipulated into a form 
allowing (direct) linear regression of the unknown param- 
eters. In this method, the single nonlinear problem is seg- 
mented into a plurality of separate problems that are exactly 
linear, thereby enabling the application of existing powerful 
linear regression algorithms such as recursive least squares. 
The individual linear sub-problems contain unknown 
parameters other than those that are identified; said param- 
eters are initially set at their nominal values, and are 
subsequently updated by the other ID processes, which are 
running concurrently. With all sub-problems sharing their 
results following each update, the results rival those of more 
computationally intensive nonlinear optimization algo- 
rithms. As long as reasonable initial parameter estimates are 
available, the approximation error for the spacecraft mass- 
property ID example is significantly smaller than that cre- 
ated by other un-modeled system parameters. 

BRIEF DESCRIPTION OF DRAWINGS 

Not Applicable. 

DETAILED DESCRIPTION OF THE 
INVENTION 

This section first describes the general aspect of the 
invention, followed by the specific spacecraft application. 
The general aspect of the invention, referred to as multiple 
concurrent recursive least squares identification (or MCRLS 
ID), is presented first using a very simple math example to 
highlight the novelty of the approach, and then in general 
form to demonstrate its generality. 

In least squares (LS) ID, unknown parameters are iden- 
tified using algorithms in which measurement data is fit to 
the underlying governing equations such that the identified 
parameter values minimize the squared error (where error is, 
for example, measurement data minus the ideal measure- 
ment data that would occur with zero noise and using the 
identified parameter values). 

As reported in the prior art, the standard form for a linear 
least squares problem, referred to as “regression form” is 
given as 

Ax=b+e (!) 

or, equivalently, 

Ax.b (2) 

where b is a vector of (perfect) measurements, e is a 
vector of measurement noise, x contains the parameters to be 
identified, and matrix A contains known variables system 
parameter values (i.e., A is noise-free). The = in the Ax=b 
representation indicates that the left and right sides of the 


Weighted, batch algorithm: x={A r WAY l A T Wb (4) 

W is a diagonal weighting matrix. Either of these algo- 
rithms can be made recursive, and the weighting matrix, W, 
15 can be chosen to weight the data according to an exponen- 
tially decaying function — as is commonly done when imple- 
mented recursively. The recursive and batch solutions are 
identical since they minimize the same cost function. 

In practice, many times the governing equations do not 
20 immediately fit exactly the form Ax=b+e, with, for example, 
noise being present in the A matrix and the x not being 
immediately and linearly separated from A and b as required. 
So the basic approach to LS ID is to find some governing 
equations (the equations of motion, for example) that con- 
25 tain the parameter values to be identified and measurement 
data. Then these equations are manipulated to conform to the 
Ax=b formulation, possibly requiring approximations along 
the way (dropping higher order terms, for example). 

In some problems, such as the spacecraft example, it is not 
30 possible to drop the nonlinear terms, since they are signifi- 
cant. A simple example that highlights this situation is 
presented by the following governing equation: 

C l X 1 +C^ 2 +c 12 x l x 2 , ‘b (5) 

35 Where b is the measurement, c 15 c 2 , and c 12 are known 
values, that vary from measurement to measurement, and x : 
and x, are the unknown parameters to be identified. This 
problem cannot be put into the form of Equation 2 , where the 
A matrix does not contain x. One approach that is feasible 
40 for a problem of this simplicity is to let 

A=[c l c 2 c l2 ];x=[x l x 2 x l x 2 ] T (6) 

With the equation now in regression form, x could be 
solved readily. However, depending on the noise present, the 
45 third element is unlikely to equal the first times the second. 
One approach would be to ignore the third ID element, but 
this is throwing away information. 

The present invention would solve this problem by 
accepting the fact that it cannot fit directly into regression 
50 form, and breaking it into two parts that can directly fit. iq 
is to be ID’ed assuming that x 2 is perfectly known, and vice 
versa. So the equations are re-written as: 

(Cl+C 12 X 2 )*l=*— C2*2 

55 (c 2 +c 12 q)x 2 =t>-c 1 ^! (7) 

or, in matrix form, 


60 

c l + C 12%2 0 

\ x l 

b — C2X2 


0 C 2 + Ci2*\ 

J 

b — ci%i 


The first equation is set up to ID x,, and treats x 2 as if it 
65 were a known quantity, substituting in the best estimate of 
x 2 , x 2 , x l5 can now be solved directly. The second equation 
does the converse. If iq and x 2 are estimated recursively, 
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with updated estimates shared with the other ID, and if the 
initial estimates are sufficiently accurate, the individual IDs 
converge correctly. It would be possible to implement this 
using a single Recursive least squares (RLS) ID using 
Equation 8, or two RLS IDs using Equation 7. The latter 5 
method is the preferred embodiment. 

One test for the error introduced by tins approximation is 
to run a simulation where in one case the secondary param- 
eters (e.g., x 2 in the top equation) are set to their true values, 
and in a second case they are set by the other ID process, and to 
to compare the accuracies of the ID results. In many real 
applications, including the spacecraft example, this approxi- 
mation error is insignificant as compared to system distur- 
bances and other unknown parameters in the system. In any 
real problem, there are always some additional unmodeled 15 
effects, and as long as the approximation error is signifi- 
cantly smaller than the error introduced by ignoring these 
effects, the approximation should be considered. 

It should be clear to those skilled in the art that this 
approach can be extended to cases where the vector of 20 
unknown parameters, x, is divided into an arbitrarily large 
number of groups, each containing an arbitrarily large 
number of parameters. The preferred embodiment is to use 
as few groups as possible to preserve the individual exact 
linear solutions to the extent possible. 25 

Another extension is that the ID algorithms used by each 
group could be something other than RLS, and could still fit 
directly into the structure presented. 

Equations describing an arbitrary number of arbitrarily 
sized groups is as follows, where the subscripts indicate the 30 
number for the group of parameters, and the groups can be 
arbitrarily sized: 

Alfa, ■■■ , *2, 111 ) X 1 = kl(*2> ■” , *n, ”■ ) 33 

■" ■")*2 =bl(X l,*3. ••• <X„. ■■■) 


A„(X 1 , ■■■ = b„(xi , ■■■ , Vi. ■■■) 40 

The first line in this equation is used as the regression 
equation in a RLS ID solution of the vector of unknown 
parameters, Xj. And so on for each group up to n. 45 

When performing RLS ID, as covered in the prior art, the 
accuracy of the ID depends on the selection of good values 
for the initial parameter estimates and the estimated cova- 
riance for the error in those estimates. It is also important to 
weight each subsequent measurement appropriately, first 50 
according to the estimated covariance of the measurement 
error, and second to exponentially weight the data so that 
newer measurements are considered more important than 
older ones (this is not strictly required, but is commonly 
done to allow the ID to track true changes in system 55 
parameters, and is a trivial adjustment to the RLS ID 
algorithm to implement). When using the present invention, 
careful selection of these parameters is even more important, 
since the IDs are dependent on each other. 

As a practical consideration, it is common for the 60 
unknown parameters in a system to have widely varying 
degrees of uncertainty. For example, in the spacecraft 
example, due to the accuracy of ground-test equipment, 
center of mass is known with better relative accuracy than 
the diagonal terms in the inertia matrix, which are known 65 
better than the off-diagonal terms. Thruster locations are 
known very well, thruster directions are less certain, and for 
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some systems, the thruster magnitude is even less certain. 
The capability of the present invention to accommodate 
these problem characteristics is an important one. 

An additional benefit of the present invention is that some 
measurements are more directly related to some parameters 
than others, and tins can be accommodated. For example, in 
the spacecraft example, when thrusters are fired to produce 
a pure torque, the resulting rotational motion as measured by 
the gyros is independent of the center of mass. Any update 
to the center of mass ID using this data is based on noise 
rather than physics and should be avoided. With the ID 
already segmented, this and other similar steps are easily 
implemented. 

If the uncertainties in the initial estimates are set too high, 
and the measurements are noisy, there is a chance that one 
of the IDs will diverge initially. The dependence of the other 
IDs on this makes it especially important. So it may be 
helpful to perform an outlier check on the measurements and 
prevent this from happening. If the estimate error covari- 
ances are set properly and the measurements are free of 
outliers, this is not a concern. 

Note that even with careful P matrix initialization and 
measurement weighting, MCRLS ID can only approach the 
mathematical optimality of a Kalman Filter or the results of 
a nonlinear optimization. However, for many applications, 
this sub-optimality is insignificant in the presence of the 
numerous other approximations and dropped terms. 

.Another benefit that is especially valuable for systems 
requiring on-line implementation on embedded processors 
where code size may be limited is that the code is very small. 
The RLS implementation consists of a small number of lines 
of code, due to the simplicity of the RLS algorithm (from 
prior art), and it can be shared by the individual IDs, leading 
to a significant savings in code size, and the time and cost 
for development, testing, and sustaining engineering. 

Detailed Description — As Applied to Spacecraft Mass-Prop- 
erty ID 

The ID method is described in detail as follows, as applied 
to spacecraft mass- and thruster-property ID. Unless other- 
wise stated, the specifics of the method refer to the preferred 
embodiment. 

The rotational and translational equations of motion for a 
thruster controlled spacecraft can be written as follows, 
where co is a 3-by-l vector containing the angular velocity 
of the body-fixed frame with respect to an inertial reference 
frame; I, is a 3-by-3 matrix containing the spacecraft inertia 
tensor (also, dyadic, matrix), measured about the true center 
of mass; L is a 3-by-n matrix containing x-y-z location of 
each thruster in the body frame; n is the number of thrusters; 
D is a 3-by-n matrix containing unit vectors indicating the 
direction of thrust in the body frame; S k is a scalar containing 
thruster magnitude scale factor applied to all thrusters. 
Includes effects of blowdown and the reduction in thrust 
when multiple thrusters are fired; F nom is a n-by-n diagonal 
matrix containing nominal strength of each thruster at full 
tank pressure; F tias , F btas , F bias (true, estimated, and error 
variants given here) is a N-by-N diagonal matrix containing 
constant off-nominal strength of each thruster at full tank 
pressure, F b ,^=t Uas +P bias i F rrmdomk is a N-by-N diagonal 
matrix containing pulse-to-pulse off-nominal strength of 
each thruster at full tank pressure; T* is a n-by- 1 vector of 1 ’ s 
and 0’s containing effective value for which thrusters fire at 
time step k, accounting for transient effects; x disturb is a 
3-by-l vector containing sum of all torques on the vehicle 
resulting from other sources (drag, gravity gradient, sepa- 
rately modeled known thruster anomalies, CMG, RWA, and 
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other calculable dynamic effects, etc.); x body is a 3-by-l 
vector containing vehicle translational acceleration, mea- 
sured in the body frame; m is a scalar containing total 
vehicle mass; and 

/■body 
J disturb 


is a 3-by-l vector containing the net force on the vehicle 
through the center of mass, with the body superscript 
indicating measurement in the body frame: 

<i>=I~ i (,(,LxD)S k (F nom +F bla ,+F ra „ d „„j [ )T k +x dl , arb - x(Im)) (10) 


tfody = m -l( DSk{Fnom + Fb . as + Frandomjc )T k + )£*) (U) 

20 

The rotational equation, contains all of the parameters we 
would like to identify, including CM location (contained in 
L), inertia (and its inverse), and Y bias . Unfortunately, the 
parameters multiply one another, and cannot be manipulated 
into the desired linear form, Axsb. One approach would be 25 
to develop a nonlinear, gradient-based optimization of these 
parameters using the full equations of motion. While pos- 
sible, that approach would be significantly more complex 
and may not produce better results than the approach taken 
here, which is to segment the identification problem into 30 
multiple sub-problems which each allow closed form solu- 
tion of the least squares problem. 

The thruster bias, CM location, inertia, and inverse- inertia 
will be identified individually by concurrently running RLS 
IDs according to the MCRLS ID approach. Note that 35 
depending on the application, not all of these values will be 
uncertain enough or critical enough to warrant ID. I -1 will 
not exactly equal the inverse of I. The equations work out to 
make I more accurate, but all IDs are presented here. 

40 

The IDs are initialized with the best prior parameter 
estimates (e.g., the nominal values). The respective estimate 
error covariance matrices are set according to the confidence 
in the initial nominal values. Then the RLS updates include 
weighting according to the sensor error covariance matrices. 45 
At each RLS update, a given ID will use the most recent 
estimates for other parameters being ID’ed (for example, 
CM ID would use the most recent estimates for inertia and 
thruster bias). The overall approach may be extended to 
other parameters in the governing equations (for example, 5Q 
the thruster directions) if they have sufficient variability and 
impact as to justify ID. 

MCRLS relies on the coupling between different param- 
eter variations to be small, so if thruster bias and inertia ID 
are tightly coupled, and they both have high variances (e.g., 55 
parameter values can be +/-50%), this approach may not 
work well. 

There is a component of the vehicle mass properties that 
can be calculated fairly well: the change in mass as fuel is 
depleted. The effect of this on the change in mass properties 60 
is calculated and is referred to as bum-time-integration 
(BTI). BTI information after each thruster firing is used to 
update the nominal mass properties. The mass ID is 
designed to ID the difference between true and nominal mass 
properties. Hie thought is that if the BTI is not perfectly 65 
accurate the mass ID will partially account for that. It should 
be clear to those skilled in the art that this approach can be 
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extended more generally to MCRLS ID where an aspect of 
one or more unknown parameters changes in a known 
manner. 

I and I -1 are symmetric, so instead of 9 free parameters in 
each, there are only 6 in each. The identifications are 
designed to directly identify these 6 parameters. 

When a pure torque is applied (i.e., the net force is zero), 
the resulting angular motion is independent of CM location, 
so the accuracy of inertia ID updates is not affected by 
inaccuracy in the estimated CM location. If the ID error 
introduced by this un-modeled coupling is a concern, inertia 
ID updates may be made only during times when a pure 
torque is applied (following the philosophy that fewer clean 
data points may be better than more noisy ones — and the fact 
that the CM-estimate error would bias the estimate). 

If using angular measurements only (e.g., with gyros), the 
CM is observable only when net translational thrust is 
applied (this is a corollary to the previous statement), so the 
CM ID is updated only at these times. 

The following algorithms enable on-line ID of inertia, 
inverse inertia, and center of mass using rotational measure- 
ments only, as would commonly be available with gyros. 
Angular acceleration is required as a minimum, and the use 
of angular rate as well allows accounting for the gyroscopic 
term (cox(Ioj)) in the governing equations. Thruster strength 
ID can be performed using rotational measurements alone, 
or with the additional use of translational acceleration mea- 
surements. 

The center of mass, C, determines the origin of the body 
frame, and thereby determines the value of L, which con- 
tains the locations of each thruster in the body frame. 
Similarly, A c , the difference between actual and nominal 
values of C determines L. A c is the value that will be 
identified here. .1 is assumed to be known perfectly; C nom is 
assumed to be known perfectly, including continual updates 
based on bum-time integration; L nom can then be calculated 
directly using these; C, and therefore L and A c also, is the 
true quantity which cannot be known perfectly; A c is the 
identified quantity, which then leads directly to C and L. 
This is done similarly for inertia and inverse inertia. 


Ac A C- C nom 

(12) 

j> 

n> 

l 

(13) 

A r iA r l -rj m 

(14) 

L = Lnom - Ac [ 1 1 ■ ■ ■ 1 ] 

(15) 


Following the MCRLS method described earlier, the 
rotational EOM is manipulated into forms that will enable 
the required IDs, individually. In the following, the values of 
the A, x, and b variables from the general regression form 
are listed for each individual ID. A few variables are defined 
now to represent quantities that appear in more than one 
regression equation. 


a k ^(Lnom X D)Sb ( + Fbias ) 2) + ^disturb 

Ck 4 DS k ( F nam 4 - Fbias)Tlc 
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-continued 

dk&ak + Ck X Lc — o>x(/o>) 


(18) 


h=dk~ Inom&k 


-continued 
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The regression form equation for center of mass ID using 
gyros only is: 


The regression form equation for inertia-matrix ID using 
gyros only, where the gyroscopic term is treated as signifi- 
cant is: 


■ 0 

~ c k, 3 

Ck, 2 


U 

<1 

Ck, 3 

0 

—Ck, 1 

, X = 

A c,2 

. ~Ck, 2 

Ck, 1 

0 


u 

<1 


bt =&-l 1 [a, - 6 j y (/&)) 


( 19 ) 


The identified value for A c is then added to the nominal 
center of mass, to give the center of mass estimate. As 
mentioned earlier, and as is present in these other IDs, ID of 
the deviation from nominal facilitates the integration of 
known changes in parameters, such as that due to the 
bum-time integration of fuel mass change and its effect on 
the nominal mass properties. 

C=C nom + A c (20) 

The regression form equation for inverse-inertia-matrix 
ID using gyros only, where the gyroscopic temi is treated as 
a disturbance is: 


A t = 



d)l 

-d) 2&3 

^>2^3 

0) 2 — d>j, 

(1)2, + (1) l (1) 2 

„2 „ 2 
^2-^3 

15 

d>ltt>3 

(1>2 

— Ct>i d>3 

&>l +- d>2^>3 

f 2 f 2 

d>3 — 

6l> 3 - d>i d>2 


—OJ\ (bo 

GJ\ (1)2 

d>3 

a.2 „2 

&1 -(i) 2 

d>l — 6l>2<^3 
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35 This is the preferred embodiment for identification of 
inertia properties since it is the only one of the three to 
directly address the gyroscopic term, and it can be updated 
during coasting rotations. However, if the spacecraft sym- 
metry and typical angular rates make the gyroscopic term 
40 negligible, there is no advantage over the other approaches. 

The regression form equation for thruster-magnitude ID 
(22) using rotational and translational measurements is: 


By treating the gyroscopic term as a disturbance, we do 
not learn anything from it; we only try to keep it from 45 
negatively impacting the ID. This means that information is 
only present when thrusters are firing, and no infomiation 
comes during coasting periods. Depending on the spacecraft 
inertia properties and typical angular rates, this temi may be 
significant enough as compared to noise and disturbances to 50 
warrant direct treatment in the ID, as is shown later. The 
magnitude of the gyroscopic term is greater for asymmetric 
spacecraft, and for high angular rates about more than one 
principal axis. 

The regression form equation for inertia-matrix ID using 55 
gyros only, where the gyroscopic term is treated as a 
disturbance is: 
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A k '■ 


b k = 


(L x D)S k diag(T k ) 
DS k diag(T k ) 

Id) 4- U) X (id) 
mx baIy 


x = diag(F bias ), 


- A k diagiF^) 


( 26 ) 


L=J-(C nom +kc)[\ 1 ... 1] (27) 

If both gyros and accelerometers are available and used 
for ID, they should be weighted accordingly (e.g., based on 
the relative accuracy in translational vs. rotational accelera- 
tions). 

The equations derived represent the relationship between 
variables at each sample update. So at each k update, and for 
each ID to be updated, the corresponding A k matrix and b k 
vector are created and run through the individual RLS ID. 

As with any system ID, the result is only as good as the 
data used. In particular, the spacecraft example is dependent 
on accurate values for angular acceleration, t%. Special care 
should be taken to estimate that as accurately as possible. 
Similarly, several aspects of the thruster model are assumed 



US 7,110,915 B2 


11 

to be perfectly known (directions, scale factor, locations) — 
they should be calibrated carefully since imprecision can 
lead to biased ID results. 

Having shown preferred embodiments for the general 
form of the invention and as applied to a specific application, 
those skilled in the art will realize many variations are 
possible which will still be within the scope and spirit of the 
claimed invention. Therefore, the scope of the invention is 
to be limited only as indicated by the following claims. 

The invention claimed is: 

1. A computer program product embodied on a computer 
readable medium having stored thereon instructions which 
when executed will perform a method for estimating a set of 
unknown parameters in a system having a set of governing 
equations describing its behavior that cannot be put into 
regression form with said unknown parameters linearly 
represented and measurements that are sampled regularly, 
the method comprising: 

(a) segmenting the set of said unknown parameters into a 
plurality of predetermined subsets wherein each indi- 
vidual subset of said unknown parameters is isolated 
linearly by manipulation of said equations; and 

(b) estimating concurrently and independently each said 
individual subset of said unknown parameters, thereby 
enabling simultaneous estimation of all unknown 
parameters, 

wherein said unknown parameters are treated as being 
perfectly known by estimating independently each said 
individual subset of unknown parameters thereby pro- 
viding unknown parameter values. 

2. The computer program product of claim 1, further 
comprising controlling when to update Unknown parameter 
subset estimates and when to ignore a set of said measure- 
ments. 

3. The computer program product of claim 1, further 
comprising enabling additional robustness to anomalous 
data by allowing nominal values, or values calculated with 
algorithms that are robust to the anomalous data, in place of 
at least one of said unknown parameter values when an 
anomaly is detected. 

4. The computer program product of claim 1, wherein said 
computer program product is implemented on-board a sys- 
tem to provide real-time parameter updates for on-board 
processes including one or more of control, estimation, 
navigation and fault tolerance. 

5. The computer program product of claim 1, wherein said 
computer program is implemented on a processor separate 
from said system, and wherein said system communicates 
measurements and results to said processor. 

6. The computer program product of claim 1, further 
comprising adjusting inputs to said system to minimize any 
negative impact on overall system performance, the method 
comprising: evaluating the unknown parameters most 
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requiring improved estimation; selecting system trajectories 
or inputs that improve the estimation for these parameters, 
while balancing other system requirements; and implement- 
ing the optimally selected system inputs. 

5 7. The computer program product of claim 1, further 

comprising calculating updated nominal values associated 
with a subset, wherein said subset comprises said parameters 
undergoing predetermined changes in time; estimating the 
deviation of said parameters from their nominal values; and 
to identifying said unknown parameters when said nominal 
values are updated. 

8. The computer program product of claim 7, wherein the 
system is a thruster-controlled spacecraft, and said changes 
are derived from the depletion of fuel mass. 

15 9. The computer program product of claim 1, wherein said 

system is a spacecraft, the unknown parameters comprising: 
the six independent parameters of the spacecraft inertia 
matrix, as measured about the center of mass; the six 
independent parameters of the inverse of said inertia matrix; 
20 and the magnitude of thrust produced by each thruster. 

10. The computer program product of claim 9, wherein 
vehicle inertia properties and motions generate both insig- 
nificant and significant gyroscopic effects, said significant 
gyroscopic effects enabling identification of said inertia 

25 matrix to continue during periods of coasting between 
thruster firings. 

11 . The computer program product of claim 9, wherein the 
spacecraft is thruster-controlled. 

12. The computer program product of claim 9, wherein 
30 the computer program is implemented on a main spacecraft 

processor, wherein said main spacecraft processor is on- 
board the spacecraft. 

13. The computer program product of claim 9, wherein 
the computer program is implemented on-board the space- 

35 craft, the spacecraft having a main processor and a second- 
ary processor communicating with the main processor, 
wherein the method is stored and run on the secondary 
processor. 

14. The computer program product of claim 9, wherein 
40 the computer program is implemented on a main spacecraft 

processor, wherein said main spacecraft processor is on- 
board the spacecraft further comprising performing at a 
ground station said calculations based on telemetry from the 
spacecraft and transmitting the calculated results back to the 
45 spacecraft. 

15. The computer program product of claim 9, wherein 
said measurements are obtained from sensors, wherein the 
sensors are at least one of gyroscopes, accelerometers, star 
trackers, sun sensors, horizon sensors, video cameras, direc- 

50 tional antennae, radar or other measurement devices that 
directly or indirectly measure spacecraft motion. 



